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(54) Procede de detection conjointe d'un ensemble de codes CDMA 



(57) L'invention concerne un procede de detection 
conjointe d'un ensemble de codes CDMA recus sur une 
pluralite d'antennes d'un recepteur. Les codes sont 
transmis par I'intermediaire d'un canal de transmission 
de matrice de transfer! A verifiant la relation e=A.d+n, d 
designant Tensemble des symboles de codes a trans- 
mettre et n un vecteur de bruit additionnel, e designant 
I'ensemble des echantillons regus. 



Le procede consiste, pour une sous-matrice de 
transfert A s correspondant a un nombre S de lignes de 
blocs de la matrice de transfert A, a determiner un en- 
semble de filtres egaliseurs permettant d'estimer sue- 
cessivement tous les symboles emis par les differents 
codes CDMA a partir d'un nombre determine d'echan- 
tillons recus e s . 

Application a la telephonie mobile. 
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Description 

[0001] L'invention concerne un precede de detection conjointe d'un ensemble de codes CDMA, utilises en telephonie 
mobile. 

5 [0002] L'UMTS, pour Universal Mobile Telecommunications System, designe la nouvelle norme de telephonie mobile 
destinee a remplacer la norme GSM. Les premieres applicatbns commerciales de cette nouvelle norme devraient voir 
le jour entre 2000 et 2005. 

[0003] Vis-a-vis de la norme GSM precitee, les principales ameliorations sont les suivantes : alors que dans le sys- 
teme GSM chaque station de base ne peut dialoguer qu'avec huit mobiles avec un debit maximum de 13 kbits/s pour 
10 chacun de ces mobiles, le systeme UMTS a pour objet, d'une part, d'augmenter le nombre d'utilisateurs, et d'autre 
part, de pouvoir allouer a chacun de ces mobiles des ressources permettant de transmettre des informations a des 
debits compris entre 4 kbits/s et 2 Mbits/s, en fonction de la nature de I'inf ormation transmise telle que parole, donnees, 
video-telephonie, hi-fi audio ou video. 

[0004] Dans ce but, deux solutions ont ete retenues au niveau du continent europeen en ce qui concerne la technique 
15 d'acces radio multi-utilisateurs : la technique dite Wideband-CDMA en mode FDD, pour Frequency Division Duplex, 
duplex a division de frequence, et la technique dite TDMA/CDMA avec detection conjointe en mode TDD pour Time 
Division Duplex, duplex a division temporelle. 

[0005] Un inconvenient majeur de la technique TDMA/CDMA en mode TDD est la grande complexity de la detection 
conjointe des codes CDMA, la detection conjointe etant definie comme le fait de proceder simultanement a la detection 
20 de tous les codes CDMA simultanement actils et recus sur plusieurs antennes, contrairement au filtre en rateau [rake 
receiver eu langage angfo-saxon) utilise habituellement en CDMA qui permet de detecter un code a la fois, les autres 
codes etant consideres comme du signal brouilleur. 

[0006] La presente invention a pour objet de remedier aux inconvenients precites par la mise en oeuvre d'un precede 
de detection conjointe de codes CDMA simplifie. 
25 [0007] En particulier, un objet de la presente invention est la mise en oeuvre d'un precede de detection conjointe de 
codes CDMA dans lequel un processus successif est mis en oeuvre, chaque appel successif permettant la mise en 
oeuvre d'une detection conjointe simplifiee d'un symbole pour chaque code CDMA simultanement actif, I'appel suc- 
cessif permettant la detection de I'ensemble des symboles de tous les codes CDMA simultanement actifs, selon une 
detection conjointe. 

30 [0008] Le procede de detection conjointe d'un ensemble de K codes CDMA simultanement actifs recus sur les Ka 
antennes d'un emetteur recepteur mobile, conforme a I'objet de la presente invention, s'applique a des codes transmis 
par I'intermediaire d'un canal de transmission lineaire multi-trajets avec bruit additif blanc gaussien verifiant la relation 
e - A.d+n ou A designe la matrice de transfert de ce canal lineaire stationnaire multitrajet, d designe I'ensemble des 
symboles avant etalement des differents codes a transmettre, n un vecteur de bruit additionnel et e I'ensemble des 

35 echantillons recus. Une estimation 5 de I'ensemble des symboles des codes est effectuee par une operation de pseudo 
inversion de la matrice de transfert de la forme b = (A + .A)" 1 .A + .e ou le symbole + designe I'operateur de transposition 
et conjugaison complexe. Le procede objet de l'invention est remarquable en ce qu'il consiste, pour une sous-matrice 
de transfert A s correspondant a un nombre S de lignes de blocs de la matrice de transfert A, a determiner un ensemble 
de filtres egaliseurs permettant d'estimer successivement tous les symboles emis par les differents codes CDMA a 

40 partir d'un nombre determine d'echantillons recus e s . 

[0009] II sera mieux compris a la lecture de la description et a I'observation des dessins ci-apres, dans lesquels, 
outre les figures 1a a 1d relatives a I'art anterieur : 

- les figures 2a et 2b represented un diagramme illustratif du procede de detection conjointe d'un ensemble de 
45 codes CDMA objet de la presente invention ; 

- la figure 3 represente, a titre illustratif, un organigramme de mise en oeuvre preferentiel du procede de detection 
conjointe d'un ensemble de codes CDMA conforme a I'objet de la presente invention. 

[0010] Prealablement a la description proprement dite du procede de detection conjointe d un ensemble de codes 
50 CDMA conforme a I'objet de la presente invention, differents elements techniques de rappel du processus de trans- 
mission de frames TDMA et de detection conjointe de ces trames seront donnes dans la description en liaison avec 
les figures 1a ; 1b, 1c et 1d. 

[001 1 ] En reference a la figure 1 a, on rappelle qu'une trame TDMA presente une duree de 1 0 ms et est decompose© 
en seize intervalles temporels, encore designes par time-slots en langage anglo-saxon. Dans chacun de ces intervalles 
55 temporels, il est possible de transmettre jusqu'a dix salves de donnees, ou bursts, grace a une separation par des 
codes CDMA orthogonaux. L'introduction des codes CDMA apporte un etalement d'un facteur 16. 
[0012] Un intervalle temporel represente 1250 temps-ch/p, la notion de ch/pen langage anglo-saxon designant un 
processus de sur-echantillonnage qui permet apres multiplication par un train binaire de realiser un etalement de 
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10 



15 



20 



25 



30 



spectre. 

[0013] Ainsi que represents en figure 1b. le processus precite autorise deux constructions de salves suivantes, 
notees Burstl et Burst2. Chaque construction comprend une sequence d'apprentissage pour estimer le canal, un 
temps de garde note G en fin de salve, et deux blocs de donnees qui representent respectivement 2x28 et 2x34 
symboles. 

[0014] On suppose que K codes CDMA sont simultanement actifs. Chaque code CDMA permet la transmission de 
2xN symboles de donnees dans une salve, confer figure 1 b. Chaque symbole de donnees de duree T s est etale par 
le code CDMA de longueur Q, encore designee par facteur d'etalement. 

[0015] Ainsi, pour un symbole de donnees, Q chips de duree T c =T S /Q sont envoyees par le canal radiofrequence. 
[0016] Le recepteur mobile est considere equipe de Ka antennes, le nombre d'antennes etant au moins egal a 2. 
[0017] Ainsi, les chips transmis sont recus au travers de Ka canaux radiofrequence, lesquels peuvent etre caracte- 
rises par le vecteur colonne de leur reponse impulsionnelle note h< ka > (ka=1 ...Ka) : cette reponse impulsionnelle etant 
echantillonnee au temps chip decrit precedemment dans la description. 

[001 8] Le tableau 1 ci-apres resume les differents parametres qui influent sur la complexity de la detection conjointe. 

TABLEAU 1 



NOM 


LIBELLE 


VALEUR 


K 


Nb de codes utilises dans rintervalle temporel considere 


1 a 10 


N 


Nb de symboles utiles dans M? burst 


28 ou 34 


Q 


Facteur d'etalement 


16 


Ka 


Nb d'antennes 


2(mobile) et + (station de base) 


W 


Longueur du canal en temps-chip 


1 a 120 


P 


Longueur du canal en temps-symbole 


1 a8 


S 


Nb d'intervalles par trame concernant I'utilisateur considere 


1 a 8 


T c 


Temps-chip 


0,462 us 


T s 


Temps-symbole (=16 T c ) 


7,38 us 


TS 


Temps-slot (intervalle) (-1250 T c ) 


577 us 


TT 


Temps-trame (=16 TS) 


10 ms 



35 



40 



[0019] On designe par W la longueur de ces vecteurs colonne precites, la longueur de ces derniers etant encore 
designee par longueur du canal en temps-chip. 

[0020] On designe par d le vecteur colonne temporellement ordonne des NxK symboles de donnees a transmettre 
dans une demie salve. Ainsi, le vecteur colonne d est obtenu en classant les donnees avant I'etalement en mettant 
en tete les premiers symboles de chaque code, puis les deuxiemes de chaque code et enfin les Niemes de chaque 
code. Les K premiers symboles du vecteur colonne d sont les premiers de chacun des K codes, ensuite les deuxiemes, 

etc. 

[0021] On designe par e le vecteur colonne des echantillons recus sur toutes les Ka antennes en ordonnant en tete 
les premiers echantillons sur les Ka antennes, puis les deuxiemes echantillons, etc. 

[0022] Apres passage dans le canal global de transmission, c'est-a-dire sequence d'etalement + canal radiofrequen- 
ce, chaque symbole de donnees emis contribue a la reception de Q+W-1 echantillons sur chaque antenne. Pour le 
kieme code CDMA, avec 1 < k < K, le transmission d'un symbole de donnees a travers le canal global precite conduit 
a la reception de Q+W-1 echantillons sur chaque antenne consideree. De cette maniere, le canal global peut etre 
considere comme un canal a Q+W-1 trajets dont la reponse impulsionnelle est le vecteur colonne verifiant la relation : 



50 



(ka.k) = h (M* c (k) 



55 



Dans la relation precedente, c< k > designe la sequence correspondant au kieme code CDMA transmis. 
[0023] Chaque symbole etant emis sur une duree de Q temps-chip, on obtient ainsi NxQ+W-1 echantillons sur cha- 
que antenne equipant le recepteur. 

[0024] La somme des contributions de quatre symboles successifs notes Symbolel, Symbole2, Symbole3 et 
Symbole4. pour N=4 symboles, W-1=2 et P=3 est representee en figure 1c. 
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[0025] En reference au formalisme introduit precedemment, il est alors possible d'ecrire I'equation matricielle ci- 
apres : 

e = Ad + n 

5 

Dans cette relation, A designe une matrice de transfert du canal lineaire stationnaire multitrajet precedemment men- 
tionne relatif au canal global, et n un vecteur colonne de bruit additif blanc Gaussien. 

[0026] Dans le but de faire apparaitre des proprietes de symetrie de la matrice A precitee, le nombre total d'echan- 
10 tillons regus doit etre un multiple de Q facteur d'etalement.Pour ce faire, des echantillons nuls peuvent etre ajoutes, 
ce qui conduit a I'ajout de coefficients nuls a la fin des b( kak ). Pour P designant la partie entiere de (W-1)/Q, P repre- 
sentant la bngueur du canal en temps-symbole, grace a introduction des echantillons nuls et des coefficients nuls 
de la matrice A considered (N+P).Q echantillons sont recus sur chaque antenne. 

[0027] La matrice A de taille [(N+P).Ka.QxN.K]est representee en figure 1d et est constitute par des sous-matrices 
is Bj de dimension [Ka.QxK] construite de la maniere ci-apres : 

- chacune des K colonnes de chaque sous-matrice B 0 a B p representee en figure Id correspond aux K a canaux 
globaux, codage + canal radiofrequence, entrelaces La notion d'entrelacement correspond a Cordonnancement 
des premiers coefficients pour chaque antenne, puis des deuxiemes, puis des troisiemes, etc. vers les Ka anten- 
na nes. En etfet chaque symbole emis apporte une contribution a Q+W-1 echantillons sur chaque antenne, soit au 
total Ka. (Q+W-1) echantillons. Ce qui, en prenant en compte quelques contributions nulles, donne Ka.Q. (P+1). 

[0028] Ainsi, la kieme colonne de chaque sous-matrice colonne relative a chaque code compris entre 1 et N est 
don nee par la relation : 



25 



30 



45 



50 



(U) ,(Ka,k) <l,k) (Ka.k) b (U> (Ka,k) 

( (b » " b » M b 2 ^2 M Q+W-1 M Q.W-, )( 0...0) T } 



[0029] Dans cette relation b^ 1 ^) . . b< Ka k > a b^^ . soulignes par une double fleche designent la contribution 

des premiers, deuxiemes et derniers elements de rang Q+W-1 a la reponse impulsionnelle du canal global multitrajets. 
[0030] La matrice de transfert A est connue apres la phase d'estimation des canaux, laquelle est effectuee a Taide 

35 des sequences d'apprentissage. La matrice de transfert A peut etre construite en calculant le produit de convolution 
de ('estimation de la reponse impulsionnelle du canal par les codes CDMA. L'objet de la detection conjointe est de 
calculer 6 une estimation de d, vecteur colonne temporellement ordonne des N.K symboles de donnees a transmettre 
dans une demi-salve, a partir de e vecteur colonne des echantillons recus sur Tensemble des Ka antennes, et bien 
entendu de la matrice de transfert A. 

40 [0031] Un processus classique de detection conjointe designe par Zero-Forcing Block Linear Equalizer pour ZF- 
BLE, propose d'effectuer le calcul de I'estimation 6 a partir de la relation : 



8 = (A + .A)* 1 .A + .e 



Dans cette relation, le symbole + en exposant represente I'operateur de trans-conjugaison, I'exposant -1 designant 
I'operateur d'inversion. 

[0032] En fait, I'estimation precitee consisle a effectuer une pseudo-inversion de la matrice de transfert A. 
[0033] Le processus d'estimation peut alors etre decompose en deux phases : 



■ Phase de desetalement permettant de se ramener aux petites tallies des vecteurs et matrices mis en jeu par le 
calcul de A + . A et de y = A + .e. 

■ Phase de resolution du systeme lineaire y =A + . A.6. La matrice A + . A etant hermitienne definie positive, la resolution 
precitee peut etre effectuee a I'aide d'une decomposition de Cholesky suivie de la resolution de deux systemes 

55 triangulares. 

[0034] La methode precedemment decrite implique une comptexite de detection conjointe critique dans la voie des- 
cendante. Le tableau ci-apres presente les resultats d'une etude de complexity dans les cas difficiles d'un reseau 
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charge et d'un canal long pour les differents services prevus pour I'UMTS : 

TABLEAU 2 



w 



15 



20 



25 



30 



35 



40 



SERVICE 


PAROLE 


384kb/s 


2 Mb/s 


Longueur canal en temps-symbole 


7 


7 


3 


Nombre d'intervalles a decoder/trame 


1 


2 


8 


% de complexite pour le desetatement 


29% 


29% 


34% 


% de complexite pour le systeme lineaire 


71 % 


71 % 


66% 


Complexite totale en MOPS 


175 


350 


570 



[0035] La complexite totale en millions d'operations, operations definies par une multiplication-accumulation com- 
plexe, par seconde, se situe a la limite des possibilites des processeurs de signaux numeriques actuels. 
[0036] Le precede, objet de la presente invention, a pour but de reduire cette complexite. 

[0037] Le procede de detection conjointe d'un ensemble de codes CDMA, objet de la presente invention, est base 
sur I'observation selon laquelle environ 70 % de la complexite totale de la detection conjointe provient de la resolution 
du systeme lineaire. En effet, I'essentiel de la complexite de la detection conjointe provient de la taille des matrices 
manipulees, taille resultant elle-meme du processus de traitement par blocs. Dans l exemple precedemment donne, 
la matrice de transfer! A est de taille 2432 x 690. 

[0038] En consequence et en reference aux figures 2a et 2b. le procede de detection conjointe d'un ensemble de K 
codes CDMA simultanement actifs recus sur les Ka antennes d'un emetteur/recepteur mobile, conforme a I'objet de 
la presente invention, consiste a considerer la transmission de ces codes par I'intermediaire d'un canal de transmission 
lineaire stationnaire multi-trajets avec bruit additionnel blanc Gaussien verifiant la relation e = A.d+n precedemment 
mentionnee dans la description, dans laquelle A designe la matrice de transfert de ce canal lineaire stationnaire multi- 
trajets ou d designe I'ensemble des symboles des differents codes a transmettre et n le vecteur de bruit additionnel, 
e designant I'ensemble des echantillons recus. 

[0039] Afin d'effectuer une estimation 5 de I'ensemble d des symboles des differents codes transmis, le procede, 
objet de la presente invention, consiste a effectuer une operation de pseudo-inversion de la matrice de transfert A, 
pseudo-inversion de la forme 8=(A+.A)" 1 .A + .e mentionnee precedemment dans la description. 
[0040] II est remarquable en ce qu'il consiste en outre, pour une sous-matrice de transfert A s correspondant a un 
nombre S de lignes de blocs de la matrice de transfert A, a determiner un ensemble de filtres egaliseurs permettan- 
td'estimer successivement tous les symboles emis par les differents codes CDMA, a partir d'un nombre determine 
d'echantillons recus. 

[0041] On indique que la sous-matrice colonne des blocs Bi avec i = 0 a p, soit la sous-matrice B 0 a B p representee 
en figure 1 d ou 2a, 2b, designe la matrice de transfert de tous les symboles emis en meme temps vers les Ka antennes. 
[0042] Ainsi que represents sur les figures 2a el 2b, pour S lignes de blocs de la matrice de transfert A, chaque 
sous-matrice de transfert A s verifie la relation (1) : 



3 s - A s d s + p 



(1) 



Dans cette relation, e s designe le sous-ensemble des echantillons recus sur les Ka antennes comprenant S blocs 
45 d'echantillons. En partant de cette relation, on peut estimer le (P+1 ) idme bloc de d s+p( ce qui correspond a un symbole 
par code CDMA simultanement actif, a condition que le nombre S verifie la relation (2) : 



(S+P).K < Ka.Q.S (2) 

Dans cette relation, on rappelle que Q designe le facteur d'etalement du code CDMA applique a chaque symbole 
transmis, P designant la longueur du canal en nombre de periodes symboles T s . 

[0043] Par decalage successif d'un bloc vis-a-vis du vecteur des symboles des codes d et du vecteur des echantillons 
e, ainsi que represents a la figure 2b, la meme sous-matrice de transfert A s permet de passer des symboles aux 
echantillons, ii est ainsi possible de construire a partir de cette sous-matrice un ensemble de filtres egaliseurs , un 
filtre par code permettant de retrouver successivement les symboles emis a partir d'un certain nombre d'echantillons 
recus sur les Ka antennes equipant le recepteur. 

[0044] Le processus de decalage represents en figures 2a et 2b peut etre utilise du debut a la fin du vecteur des 
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echantillons en appliquant une fenetre glissante d'un pas correspondant Ka.Q pour obtenir successivement tous les 
symboles emis par le code precite. 

[0045] Ainsi, en reference aux figures 2a et 2b, la meme sous-matrice A s permet de calculer successivement tous 
les blocs de la matrice des symboles pour realiser Pestimation consideree. 
s [0046] Sur la figure 2a et la figure 2b, on Indique que dgj p designe le vecteur (ligne) transpose du vecteur (colonne) 

[0047] Le procede, objet de la presente invention, reste en outre valable pour les extremites du vecteur d des sym- 
boles des differents codes. En effet, I'adjonction aux deux extremites de blocs nuls au vecteur d precite n'entraine pas 
de consequence prejudiciable, la relation precedente liant chaque sous-matrice de transfert A s , le vecteur des sym- 
70 boles des differents codes correspondant et le vecteur des echantillons e s restant valable. 

[0048] La taille minimum de la sous-matrice a prendre en compte est donnee par la relation (2) precedemment 
mentionnee dans la description et permettant d'estimer un bloc du vecteur d s+P . Cette relation indique que le vecteur 
e s doit etre plus long que le vecteur d s+P . 

[0049] Toutefois, les performances peuvent etre ameliorees en utilisant tous les echantillons auxquels le bloc a 
15 estimer du vecteur d a contribue. Ce mode operatoire permet de dresser le tableau 3 ci-apres qui represente en f onction 
de P, la plage dans laquelle il faut choisir la valeur S, la meilleure valeur etant celle qui correspond au compromis 
complexite/performance recherche. 



TABLEAU 3 



p 


1 


4 


7 


S MIN 


1 


2 


4 


S MAX 


2 


5 


8 



[0050] II taut noter que le choix de S = S MAX implique toutetois que la sous-matrice de transfert A s presente toutefois 
dix fois moins de lignes et cinq fois moins de colonnes que la matrice de transfert A. 

[0051] En outre, selon un aspect particulierement remarquable du precede, objet de la presente invention, on indique 
que pour une matrice de transfert A de dimension [Ka.Q.(N+P)xK.N] qui peut etre consideree comme une matrice de 
taille [(N+P)xN] composee de blocs de taille [Ka.QxK], les filtres muttidimensionnels egaliseurs sont formes par la 
(P+1)* me ligne de blocs de la matrice As' 1 pseudo-inverse de la sous-matrice de transfert A s et verifiant la relation : 

(A s ) "=(^sYK 



[0052] Si en outre un seul code est alloue au recepteur mobile, il suffit alors de calculer la ligne (de la ligne de bloc 
calculee A3" 1 ) correspondant au code recherche. 

[0053] Ainsi, la ligne de rang k de cette ligne de blocs definit le filtre egaliseur pour le code de rang k. 
[0054] Dans le cas ou le calcul direct de la matrice pseudo-inverse de la sous-matrice de transfert A s s'avere trap 
complexe, alors, le meme type de resolution que dans le cas de la detection conjointe classique peut etre utilise. 
[0055] D'une maniere generale, on indique que les donnees decodees sont le resultat du filtrage du produit A s + .e s 
par les lignes de (As+.As)' 1 calculees. Ainsi, (S+P).K multiplications sont necessaires pour chaque symbole de don- 
nees. Dans le cas ou un seul code est alloue au mobile, la complexity globale de cette etape est ramenee a (S+P).K. 
N multiplications, ce qui represente une quantite negligeable par rapport au processus de detection conjointe classique. 
[0056] Une description plus detaillee d'un mode de mise en oeuvre preferentiel du procede de detection conjointe, 
objet de la presente invention, sera maintenant donnee en liaison avec la figure 3. 

[0057] D'une maniere generale, on indique que le procede, objet de la presente invention, comprend une phase 
d'organisation des donnees suivie d'une phase de calcul proprement dite, la phase de calcul permettant de definir le 
processus de filtrage applique a Pensemble des echantillons recus constituant le vecteur colonne des echantillons 
recus e. 

[0058] Ainsi que represente sur la figure 3, la phase d'organisation des donnees peut comporter une etape 1000 
consistant a organiser temporellement le vecteur colonne e des echantillons recus sur les Ka antennes. Cette etape 
1000 peut etre suivie d'une etape 1001 consistant a construire les blocs Bj constitutifs de la matrice de transfert A 
precedemment mentionnee dans la description. Contormement a un aspect specifique du procede, objet de la presente 
invention, I'etape 1001 peut consister a etablir des blocs B { de facon que la k ldnie colonne soit formee par les premiers 
echantillons des reponses impulsionnelles des canaux globaux du code de rang k vers les Ka antennes, puis les 
deuxiemes echantillons et ainsi de suite, ainsi que mentionne precedemment dans la description. 
[0059] La phase d'organisation des donnees precitee, constitute par les etapes 1000 et 1001 , peut alors etre suivie 
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d'une phase de filtrage proprement dite. 

[0060] La phase de fittrage proprement dite peut comporter, ainsi que represents de meme en figure 3, une etape 
1003 consistant a choisir le nombre S de lignes de la matrice de transfer! A pour en tait constituer !a sous-matrice A s 
ainsi que mentionne precedemment dans la description. On rappelle que le choix de S peut etre effectue conformement 

5 aux criteres indiques precedemment dans la description en liaison avec le tableau 3. L'etape 1 003 comporte egalement 
une etape de critere de test du nombre d'elements de la sous-matrice A s ainsi constitute a une valeur de seuil deter- 
minee, notee N 0 . On comprend en particulier que la valeur de seuil N 0 peut etre definie de maniere arbitraire, experi- 
mentalemenL en fonction des possibilites de calcul dont on dispose. L'etape 1003 permet ainsi, suite au choix de la 
taille de la sous-matrice A s precitee et si cette sous-matrice est suffisamment petite, le nombre d'elements de cette 

10 derniere N AS etant inferieura N 0 , d'effectuer une etape 1004 consistant a calculer directement la P+1® me ligne de blocs 
de la matrice pseudo-inverse de la sous-matrice A s , cette etape 1004 etant elle-meme suivie d'une etape 1005 con- 
sistant a filtrer les echantillons recus du vecteur colonne d'echantillons e par I'intermediaire de la ligne de rang P+1 
de blocs de la matrice pseudo-inverse de la sous-matrice A s . On indique que le filtrage realise a l'etape 1005 est en 
fait applique successivement au bloc d'echantillons e s de taille Ka.Q.S par decalage successif de Ka.Q echantillons, 

15 ce decalage correspondant au mode operatoire represents precedemment dans la description en liaison avec les 
figures 2a et 2b. 

[0061] Dans le cas ou le nombre d'elements N AS de la sous-matrice A s n'est pas inferieur a la valeur de seuil 1%, 
l'etape 1003 est au contraire suivie d'une etape 1006 consistant a calculer le produit de la sous-matrice A s transcon- 
juguee, notee A s + et de la sous-matrice A s precitee. Cette etape 1 006 est elle-meme suivie d'une etape 1 007 consistant 

20 a calculer au moins une ligne de la P+1 6me ligne de blocs de I'inverse du produit A S + .A S calcule a l'etape 1006 prece- 
dente. On indique que le calcul de l'etape 1 007 peut etre realise par la mise en oeuvre d'une decomposition de Cholesky. 
[0062] L'etape 1 007 est elle-meme suivie d'une etape 1008 consistant a calculer, pour tout bloc d'echantillons recus 
es de taille Ka.Q.S du vecteur colonne d'echantillons e le produit As + .e s de la sous-matrice transconjuguee A s + de la 
sous-matrice A s et de ce bloc d'echantillons recus. De maniere semblable a l'etape 1005, on rappelle que chaque bloc 

25 e s est obtenu par decalage successif de Ka.Q echantillons ainsi qu'illustre en figures 2a et 2b precedemment dans la 
description. 

[0063] L'etape 1008 est alors suivie d'une etape 1009 consistant a filtrer chaque produit As + .e s obtenu a l'etape 
1008 precedente par I'intermediaire de la ou des lignes de la P+ieme ligne de blocs de I'inverse du produit A S +.A S . 

30 

Revendications 

1. Procede de detection conjointe d'un ensemble de K codes CDMA simultanement actifs recus sur les Ka antennes 
d'un emetteur/recepteur mobile, ces codes etant transmis par I'intermediaire d'un canal de transmission lineaire 

35 stationnaire multitrajets avec bruit additionnel blanc Gaussien et verifiant la relation e = A d + n ou A designe la 

matrice de transfer! de ce canal lineaire stationnaire multitrajet, d designe I'ensemble des symboles des differents • 
codes a transmettre et n le vecteur de bruit additionnel, e designant I'ensemble des echantillons recus, procede 
dans lequel une estimation 6 de I'ensemble d des symboles des differents codes est effectuee par une operation 
de pseudo inversion de ladite matrice de transfer! de la forme 6= (A + .A)" 1 .A + .eou le symbole + designe I'operateur 

40 de transposition e! de conjugaison complexe, caracterise en ce que ledit procede consiste, pour un ensemble de 

sous-matrices de transfert As correspondant a un nombre S de lignes de blocs de ladite matrice de transfer! A, a 
determiner un ensemble de filtres egaliseurs permettant a partir d'un nombre determine d'echantillons recus es 
d'estimer successivement tous les symboles emis par les differents codes CDMA. 

45 2. Procede selon la revendication 1 , caracterise en ce que pour toute sous-matrice de transfer! A s verifiant la relation : 

e s = A s .d s+P 

50 e s designant le sous-ensemble des echantillons recus sur les Ka antennes comprenant S blocs d'echantillons, le 

nombre S verifie la relation : 

(S+P).K < Ka.Q.S 

55 

relation dans laquelle Q designe le facteuV d'etalement du code CDMA, applique a chaque symbole transmis, P 
designant la longueur du canal en nombre de pertodes symboles. 
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3. Procede selon Tune des revendications 1 ou 2, caracterise en ce que pour une matrice de transfert A, de dimension 
[MxN KJavec M = Ka.Q.(N+P), N designant le nombre de symboles transmis par code CDMA par demi-salve 
d'emission, cette matrice de transfert etant formee de [(N+P)xN] blocs, chaque bloc etant de taille[Ka.QxK], les 
filtres multidimensionnels egaliseurs sont formes par la (P+1) e ligne de blocs de la matrice (Ag)" 1 pseudo-inverse 

5 de ladite sous-matrice de transfert A s et verifiant la relation : 

<A S )' 1 =(a^-A s )" 1 .A* 

10 le symbole -1 designant I'operateur de pseudo-inversion, la ligne de rang k de cette ligne de blocs definissant le 

filtre pour la code de rang k. 

4. Procede selon la revendication 3, caracterise en ce que pour chaque code une fenetre glissante d'un pas Ka.Q 
est appliquee sur les echantillons pour obtenir successivement tous les symboles emis par ce code. 

15 

5. Procede selon I'une des revendications 1 a 4, caracterise en ce que celui-ci comporte les etapes consistant a : 

- organiser temporeilement le vecteur colonne e des echantillons recus sur les Ka antennes ; 

- construire les blocs Bj, la kieme colonne de la matrice constitute des Bq, B v ...B p places les uns au-dessus 
20 des autres en colonne etant constitute des premiers echantillons des reponses impulsionnelles des canaux 

globaux du code k vers les Ka antennes, puis les deuxiemes et ainsi de suite ; 

choisir le nombre S de lignes et construire la sous-matrice A s a partir des blocs Bj correspondents obtenus a 
I'etape precedente ; 

- suite a un critere de choix sur la taille de la sous-matrice A s , si le nombre d'elements de la sous-matrice A s 
25 est inferieur a une valeur de seuil dtterminee, 

calculer la (P+1 ) 6me ligne de blocs de la matrice pseudo inverse de la sous-matrice A s ; 

filtrer les echantillons recus du vecteur colonne d'echantillons e par I'intermediaire de la ligne de rang P+1 de 

blocs de la matrice pseudo inverse de la sous-matrice A s ; 

appliquer successivement le filtrage realise a Tetape precedente aux blocs d'echantillons e s de taille Ka.Q.S 
30 par decalage successif de Ka.Q echantillons ; sinon, le nombre d'elements de la sous-matrice A s n'etant pas 

inferieur a ladite valeur de seuil determinee, 

a calculer le produit de la sous-matrice A s transconjuguee et de la sous-matrice A s ; 

calculer au moins une ligne de la (P+1 ) bme ligne de blocs de ('inverse du produit calcule a Tetape precedente ; 

- calculer, pour tout bloc d'echantillons recus es de taille Ka.Q.S du vecteur colonne d'echantillons e, le produit 
35 A s + . e s de la sous-matrice transconjuguee A s + de la sous-matrice A s et de ce bloc d'echantillons re$us, chaque 

bloc e s ttant obtenu par decalage successif de Ka.Q echantillons ; 

- a filtrer chaque produit As + .e s obtenu a i'etape precedente par I'intermediaire de ladite au moins une ligne de 
la (P+1)« me ligne de blocs de I'inverse du produit As + .A S . 
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